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Abstract 

A method is presented for first-principles calculations of inelastic mean free paths and stopping 
powers in condensed matter over a broad energy range. The method is based on ab initio calcu- 
lations of the dielectric function in the long wavelength limit using a real-space Green's function 
formalism, together with extensions to finite momentum transfer. From these results we obtain 
the loss function and related quantities such as optical-oscillator strengths and mean excitation 
energies. From a many-pole representation of the dielectric function we then obtain the electron 
self-energy and inelastic mean free paths (IMFP). Finally using our calculated dielectric function 
and the optical-data model of Fernandez- Varea et al, we obtain collision stopping powers (CSP) 
and penetration ranges. The results are consistent with semi-empirical approaches and with ex- 
periment. 

PACS numbers: 77.22. Ch; 34.80.-i; 34.50.Bw 



1 



I. I. INTRODUCTION 



The effect of inelastic losses on fast electrons has long been of theoretical and exper- 
imental interest^ 2 ^, and continues to be an area of active development^. Theoretical 
calculations of such losses depend on the dielectric response of a material over a broad 
spectrum. Moreover, calculations of losses at low energies are particularly sensitive to the 
excitation spectrum of a material. While first-principles approaches have been developed 
for calculations of losses at low energies, i.e., up to a few tens of e"V^, these methods are 
computationally intensive and may be difficult to implement. Thus detailed calculations 
of inelastic losses have generally been limited to semi-empirical approache o 4 ! 9 ' 10 , based on 
experimental optical data. On the other hand, experimental data over a sufficiently broad 
spectrum are not readily available for many materials of interest. 

In an effort to overcome these difficulties, we present here a first-principles, real-space ap- 
proach for general calculations of inelastic losses. The approach is applicable to both periodic 
and aperiodic condensed matter systems throughout the periodic table. Our calculations 
are based on ab initio calculations of the complex dielectric function e(oj) = e 1 (c ( j) + ^(cu) 
in the long- wavelength limit, together with extensions to finite momentum transfer—. The 
calculations of e(u) are carried out using an all electron, real-space Green's function (RSGF) 
formalism as implemented in a generalization of the FEFF8 cod o 12 ' 13 for full-spectrum cal- 
culations of optical constants. 

We focus in this paper on several physical quantities which characterize the inelastic 
interactions of a fast probe electron, a photo-electron, or other charged particle in condensed 
matter. These include the inelastic mean- free-path (IMFP) and the collision stopping-power 
(CSP). Each of these quantities depends on the complex dielectric function e(oj) through the 
loss function for a given material — [Ime _1 (c<j)] = € 2 (uj) / [ei(u)) 2 + e 2 (^) 2 ], which is calculated 
here up to x-ray energies. The loss function is directly related to the optical oscillator 
strength (OOS). From the OOS we obtain values of the mean excitation energy / which 
characterizes the distribution of excitations (e.g., plasmons, particle- hole excitations, etc.). 
Recently a comprehensive relativistic treatment of inelastic losses and scattering within 
the first Born approximation has been developed by Fernandez- Varea et alA Their semi- 
empirical approach requires experimental optical data as input, and is referred to here as the 
optical data model (ODM). This approach has the advantage that calculations of quantities 



2 



such as the CSP are reduced to a single quadrature. To facilitate precise comparisons, we 
have used their formulation for our CSP calculations, except for the substitution of our ab 
initio dielectric function. Our approach is therefore referred to as the u ab initio data model" 
(ADM). We have also compared IMFPs calculated using both the ADM and a one-particle 
self- energy approach. 

Formally the IMFP and the CSP are related to energy moments of the differential cross- 
section (DCS) for inelastic collisions da/duj of a fast probe electron (or other charged particle) 
of initial kinetic energy E with energy loss u. The inverse IMFP is proportional to the zeroth 
moment of the DCS 



where n a = N/V is the atomic number density. Here and elsewhere in this paper we use 



energies in Hartrees (H « 27.2 eV), unless otherwise specified. The CSP, here denoted by 
S(E), is proportional to the first moment of the DCS 



Since S(E) = —dE/dx, this quantity has units of force. From an integral of 1/S(E) over 
energy we then obtain the net penetration range or path length R(E). Implicit in Eqs. (JT|) 
and (J2J are the kinematics of the colliding particles. In this paper we choose kinematics 
relevant for electrons, but we could alternatively use similar equations to describe protons 
or other ions by suitably modifing the domain of integration that defines da/du (see below). 
Regardless of the probe, the sample is characterized by the dielectric function e(q, to). In 
this paper we consider cubic materials which we approximate as isotropic, i.e., in which the 
dielectric function depends only on the magnitude of the momentum transfer q — |q|. 

The DCS may be considered as the sum of longitudinal (instantaneous Coulomb) and 
transverse (virtual photon) contributions, denoted below with subscripts L and T respec- 
tively. The detailed relativistic form of the relationship between each contribution to the 
DCS and the loss function is obtained by integrating the double differential cross-section 
(DDCS) over the kinematically allowed values of momentum-transfer 4 , 




(1) 



Hartree atomic units (m = % = e 2 = 1). Thus distances are in Bohrs (a ~ 0.529 A) and 






(3) 
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where 

da(q,uj) _ da L (q,u) da T (q,uj) ^ 
dqdu dqduj dqdu 

As an example of how the dielectric function determines the DDCS, we recall the familiar 
non-relativistic result: 

da(q,u) da L (q,u) 1 _ x 

~~r~A — = — n = ~o 2 lme 5 

dqdu dqau 27m a qv z 

The relativistic analog of Eq. (JHJ) is similar and is given explicitly in Eqs. (8) and (9) of 
RefA 

One of the main goals of this work is to calculate mean excitation energies I and IMFPs 
for general condensed matter systems over an energy range up to about 100 keV. Another 
goal is to calculate CSPs and penetration ranges over a range of order 10 MeV. We compare 
our results both with other semi-empirical approaches and with experimental data and 
tabulations. 



II. II. MODEL DIELECTRIC FUNCTION 



Both the IMFP and the CSP can be computed as convolutions of the momentum-transfer 
and energy- loss dependent inverse dielectric function e~ 1 (q,u), with relativistic weighting 
functions. The precise details of the weighting functions are discussed further below. In this 
Section we discuss the extension of our ab initio calculation of e(q, to) in the long wavelength 
(q — ► 0) limit to finite g 4 i 9 i 14 i 15 i 16 . i n this work e(u) = e(0, oS) is calculated from the UV 
to x-ray energies using the ab initio real-space Green's function code FEFF80 P 1 ? 11 ? , which 
sums the contributions to the spectra over all occupied core and semi-core initial states. As 
an example, the calculated loss function for Ag is shown in Fig. ^ 

We have chosen to discuss the extension to fmite-g in terms of the loss function, but we 
could just as well have used the OOS g(ui), with differs by a factor proportional to u, i.e., 

g( U ) = -~ U 1m[e-\u)], (6) 

where = 47m a Z is the all-electron plasma frequency. As an illustration of the quantitative 
agreement of our approach, three ab initio OOS calculations, spanning a range of atomic 
numbers, are compared to experiment in Fig. El Clearly the approximations in FEFF8 
such as the use of atomic core initial states in the OOS calculation and muffin-tin scattering 
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FIG. 1: Loss function -[ImT^w)] of fee silver as calculated in this work (solid line) and from 
experimen t 18 ! 19 (dots). 

TABLE I: Mean excitation energies for several elements as calculated in this work, and for com- 
parison, results calculated from experimenta l 19 ! 21 optical constants, and recommended (ICRU^ 2 -) 



values. 

Element I theory (eV) I expt (eV) ICRU (eV) 

Aluminum 165 167^-, 1631 9 . 166 

Silicon 174 1732i 173 

Copper 312 3191 9 - 322 

Silver 420 3821 9 . 470 

Gold 662 7531 9 - 790 



potentials are adequate to yield good agreement with experiment for UV energies and above. 
Additional examples are tabulated on the WWW- . For optical frequencies and below, 
however, the agreement is only semi-quantitative, but the errors tend to be suppressed in 
the OOS due to the overall factor of u in Eq. (JHJ). 

A global measure of the excitation spectra is given by the "mean excitation energy" 
In/ = (In a;), where the "mean" (. . .) refers to an average with respect to the OOS weighting 
function, i.e., 

J dujg(uj) 

The quantity / appears in expressions for the collision stopping power, as shown in Sec. V. 
below. In Table U theoretical values of /, as calculated from our OOS spectra, are compared 
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FIG. 2: Optical oscillator strengths for a) Cu (upper), b) Ag (middle) and c) Au (lower) as 
calculated in this work (solid) using the FEFF80P code and compared to experimen t 18 ! 19 (points). 

with those calculated from experimental optical constants^ 9 , and also with internationally 
recommended (ICRU) values for several elements. For low Z elements, the theoretical values 
of / are clearly in good agreement with measured values. For high Z elements I is predicted 
by Thomas-Fermi models to be proportional to Z. The proportionality constant can be 
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determined experimentally to give a semi-empirical "rule of thumb" I ~ 10 Z (eV). In 
the high Z regime the agreement between theory and experiment appears to be only semi- 
quantitative, but it should be mentioned that the ratio of I to Z is not in fact constant but 
rather it varies from around 9 to around 11 and has been measured to be as low as 8.32 
for lanthanum^. Furthermore, only the logarithm of / is needed for the determination of 
physical quantities, and typical errors in In / from the values in in Table I are only a few 
percent. For example, for the case of 100 keV electrons in gold, the error in stopping powers 
calculated using experimental versus theoretical values of / is only around 2%. 

In the IMFP calculations presented here we consider two different extensions to finite q, 
as described below. Since our calculations show that both lead to similar results for the 
IMFP, we only present calculations of the CSP with one of these extensions. However, all 
our calculations use the same full-spectrum calculations of e{oj). For our IMFP calculations, 
we have used the many-pole representation of the dielectric function of Ref.— , i.e. we 
approximate our calculated loss function — Im [e _1 (o;)] as a sum of many (typically of order 
100) discrete poles. This "many-pole" model, denoted by e^ v 1 (c<j), has the standard analytic 
form for dielectric response, 



where 5uj is a small damping term, comparable to the pole separations. Fig. |3] compares 
the IMFP for Cu as calculated using both our many-pole model and a single plasmon-pole 
model 24 . This single-pole model is essentially an Einstein-model for the response in which 
excitations (for a given momentum transfer q) occur at the plasmon excitation energy u g . 
Thus the single-pole model is a special case of the many-pole model in which all but one of 
the weights Qj appearing in Eq. (jHJ) are set to zero. In our many-pole representation, the 
parameters ujj are taken to be evenly spaced along the energy-loss axis, and the weights gj 
are fixed by matching our calculation of Im[e -1 (u;)] according to the formula 



where the integration region Aj is from (ujj + u>j^i)/2 to (ujj + Uj + {)/2 and the similarity 
with Eq. (jUJ) is apparent. Finally, the extension to finite q is obtained by shifting the pole 
locations via the substitution^ 4 - 






(9) 




4 ' 



(10) 
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where vf = hp/ mis the Fermi velocity as calculated at the mean interstitial electron density 
from the FEFF8 code. Further details of our approach, though not essential to our discussion 
here, are given in Refs.— ,— and^. The above substitution is sufficient to induce the so-called 
"Bethe ridge" for large momentum-transfer where the loss function is peaked about the point 
u = q 2 /2. In other words, our model for large q satisfies the approximate relation 



where Q = q 2 /2. Consequently the above extension can be regarded as an interpolation 
formula between small and large Q. 

For CSP calculations, our aim here is to replace experimental optical data (which is 
used as input in the ODM of Refj») with theoretical "optical data" from our ab initio 
calculation of e(u), i.e., with an ab initio data model (ADM). Thus for consistency we follow 
the formulation of Ref.— as closely as possible in comparisons with their CSP results, In 
particular we have also implemented their delta-oscillato r 141 ^ extension to finite q for our 
CSP calculations. For non-relativistic probe electrons the delta-oscillator model extends 
e(u) to finite q according to the relation 



where Z(Q) is the number of electrons that contribute to the zero momentum-transfer sum- 
rule with upper energy limit Q, 



Because Z(Q) approaches Z for large q we see that the extension to finite q in Eq. (JT2)l gives 
the Bethe ridge in much the same way as that of Eq. fTTf . Although the finite-g extension 
algorithms in this paper differ somewhat, we do not expect our non-relativistic results to 
depend significantly on the details. The reason is that both algorithms reduce to the correct 
long-wavelength limit for low q, and both give the correct Bethe ridge dispersion for high q. 
This expectation is supported by the IMFP results in Sec. IV. Moreover our results for the q 
dependence are roughly consistent with the explicit real space calculations of — Im[e _1 (g, u)] 
at finite q of Soininen et al.— . 




(11) 




Im[e-V)]0(u;-Q) , 



(12) 




(13) 
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III. III. ELECTRON SELF-ENERGY 



Inelastic losses in the propagation of a fast charged particle can be expressed in terms 
of one-particle self-energy S(i?). This complex- valued quantity is a dynamically screened 
exchange-correlation contribution to the quasi-particle energy-momentum relation 



where p is the quasi-particle momentum. Our approach for calculating T,(E) is based on 
the U GW" approximation of Hedin 28 , together with our many-pole representation of the 
dielectric function, as summarized above^. In the GW method the vertex-corrections to the 
electron self-energy are neglected, yielding an expression for T,(E) in terms of the electron 
propagator G and the screened Coulomb potential acting on an electron W, i.e., 



where r] is a positive infinitesimal and spatial indices (x, x') have been suppressed for clarity. 
Within the RSGF approach, the propagator G is calculated using a multiple-scattering 
expansion G = Gq + GqIGq + • • •— . However, for simplicity in this work, we neglect the 
multiple-scattering terms (which would give rise to fluctuations in the self-energy) and simply 
use the free propagator Go for a homogeneous electron gas at the mean interstitial density. 
Then the screened Coulomb interaction for a spatially homogeneous system can be obtained 
from the Fourier transform W(q, u) = T [W(x — x', t)} and can be expressed in terms of the 
Coulomb potential v q = An/q 2 and the dielectric function e(q, to) as: 



The calculations of £(-E) are then carried out using the many-pole representation of Eq. (jHJ) 
and (fTT)|) . With this homogeneous model, our calculated ^{E) is then the average self-energy 
in the system. Further details are given in Ref.— . 

IV. IV. INELASTIC MEAN FREE PATH 

We first calculate Eq. for the IMFP in terms of the excited state self-energy T,(E) of 
the fast electron. 



(14) 




(15) 



W(q,u) = e 1 (q,u)v g . 



(16) 




(17) 
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FIG. 3: Inelastic mean free paths for copper calculated using the same ab initio dielectric function 
as the basis of two different theoretical models: The many-pole self-energy (MPM) model of Eq. Q17|) 
and the single-pole self-energy model (described in the text). These theoretical results are compared 
to: Exp. (aft 9 (squares), Exps. (b) (circles, the references for Exps. (b) are given in Reb 5 -), and a 
semi-emperical curve which is described in Reb 5 -. 

Eq. (|17|) for the IMFP is consistent with the decay of a single electron wavefunction whose 
time dependence is given by e~ tE<J> ^ 1 . Eq. (fTTj) is also equivalent to Eq. ((H), because the 
self-energy is proportional to the forward scattering amplitude, i.e., 

ImS(p) = -27rn a Im/(p,p) , (18) 

i.e., the equivalence of Eq. (0) and Eq. (JT7jl follows from the optical theorem. 

The explicit dependence of the self-energy on the dielectric function in Eq. (|16jh the 
many-pole model of Eq. (jHJ), and (jlUj) . and the full-spectrum FEFF80P code are all that are 
needed to carry out ab initio calculations of IMFPs according to Eq. (|17|) . Our many-pole 
model IMFP calculations (labeled MPM) are shown for several materials in Figs. El and 0J 
together with best fits 5 to currently available data. The fit lines in Figs. H3 and 0] are based on 
multiple data sets which were taken up to 3000 eV and are expected to accurately describe 
the IMFP as low as 50 eV. Fig. 0] also shows a calculation (labeled "ADM" ) which uses our 
ab initio e(u) as input data to the semi-empirical optical-data model of Ref A Note that 
the MPM and ADM results are in in good agreement with each other, which verifies that 
the different extensions to finite q discussed in Sec. II. do not lead to significantly different 
results. Both theoretical models are plotted here over the expected range of validity of the 
fit line, but can be extended with our codes to energies up to about 100 keV. Although 
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FIG. 4: Inelastic mean free paths for a) silver (upper) and b) gold (lower) calculated using the 
same ab initio dielectric function as the basis for two different theoretical models: the many-pole 
self-energy (MPM) model of Eq. (|17|). and the "a& initio data" model (ADM) described in the 
introduction. The theoretical results are compared to a semi-emperical curved and to multiple 
experimental data sets. The references for the Exps. are given in Ref£. The theoretical models 
are plotted over the expected range of validity of the semi-emperical curve. 

the agreement with experiment is reasonable, our calculations tend to underestimate the 
experimental IMFP somewhat for high Z materials. 

V. V. STOPPING POWER 

As noted in the introduction the CSP is the net reaction force S(E) = —dE/dx due to 
electronic collisions at a given energy E that slows a fast probe electron. Over the range 
of energies from about 10 eV up to about 10 MeV the CSP is the main contribution to the 





11 



total stopping power. Above this energy the total stopping power may be dominated by 
bremsstrahlung. The CSP is calculated in this work using Eq. (J2J), where the DCS is related 
to our ab initio loss function using the formulation of RefA This model is thought to be 
applicable with confidence for energies above about 100 eV, and appears to be applicable 
as low as about 10 eV. However, it is not obvious why a model based on the first Born 
approximation should be valid at such low energies. In the relativistic limit Eq. (J2J) reduces 
to the well-known Bethe-formul a 1 ' 2 ! 3 ' 32 for the stopping power 



Also appearing in Eq. (|19|) are the "mean excitation energy" / defined in Sec. II., and Fano's 
density correction 3 5p. The density correction 5p is due solely to transverse interactions, 
and can be neglected for non-relativistic particles. A detailed description of how Sp can be 
calculated as a functional of the loss function, is given in Ref.—. Fig. 03 shows the density 
correction 5f{E) for copper as calculated using our ab initio dielectric function, and for 
comparison, the semi-empirical values used by ESTAR- 3 ^. ESTAR is an on-line implemen- 
tation of Eq. (|T9"|) which semi-empirical values of / as input. The mean excitation energy 
and the density correction have, heretofore, been difficult to calculate from first-principles, 
as they require accurate values of the OOS over a very large energy spectrum. However, our 
full-spectrum approach clearly gives reasonable agreement with experiment. 

For relativistic probe particles, the excellent agreement of the Bethe formula in Eq. ()19j) 
for the CSP is well known, so we have included data from ESTAR in lieu of experiment where 
necessary. The difference between Eq. ()19|) and Eq. (J2J) only appears in the non-relativistic 
regime, and can be seen in Fig. El where Eq. (fTT?|) begins to fail around 5000 eV. In order 
to calculate CSPs that are in good agreement with experiment for both non-relativistic and 
relativistic probe electrons we apply Eq. (j2J) with the more general form of da/duj given in 
Ref.—, but using our calculated dielectric function as input. 

Figs. |H1 and [7| show that the calculated CSPs for high energy electrons for the systems 
studied here (Cu, Ag and Au) are in good agreement with the ESTAR results. For lower 



(19) 



where 7 = (1 — v 2 /c 2 ) x l 2 is the relativistic dilation factor, and ^(7) is given by 




(20) 
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FIG. 5: Fano's density effect correction to the stopping power from Eq. (|19|) as calculated in this 
work (solid), and compared to semi-empirical values^ 3 - for copper (dashes). 
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FIG. 6: Collision stopping power for copper as calculated using the ab initio dielectric function of 
this work (dashes) in the ADM (see text). Also shown are semi-empirical values of the CSP from 
ESTAR 33 (open squares), and semi-empirical CSP values (labeled TPP-calc) based on the Penn 
model 30 (circles), and CSP values from experimental data: (+) 34 , (x)^, and (solid squares)^. 

energies the CSP calculations of this work show significantly better agreement with data 
than the Bethe formula. 

It is interesting to note that the CSP falls off approximately as a power-law of the energy 
over a few decades beyond the peak loss, but before the transverse effects take over. Because 
of this, the range R(E) defined by 

R(E) =£my (21) 

is also well approximated as a power-law. This approximate power-law dependence of R(E) 
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FIG. 7: Collision stopping powers for a) Ag (upper) and b) Au (lower), with labels as in Fig. 
Also shown for Au are the semi-empirical CSP values as calculated in Ref?A (labelled FSDL-calc), 
and CSP values from experiment 34 . 

is shown in Fig. |Hl for silver, and compared with our full calculation and ESTAR. 



VI. VI. CONCLUSIONS 

We have presented a general real-space Green's function approach for ab initio calculations 
of inelastic losses and stopping powers in condensed matter. Unlike most current approaches, 
our method is based on ab initio calculations of dielectric response, and does not rely on 
empirical optical data. We find that our many-pole self-energy model, which is derived from 
our many-pole dielectric function, yields inelastic mean free paths in better agreement with 
experimental data than the single-pole plasmon model. Thus our many-pole self-energy more 
accurately accounts for inelastic losses in various electron and x-ray spectroscopies than the 
single plasmon-pole self-energy. We also find that accurate calculations of inelastic losses 
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FIG. 8: The range R(E) as given by Eq. (|21j) of electrons in silver as calculated in this work (solid) 
and compared to ESTAR (squares) and to a pure power- law R(E) = 0.271E 5 / 3 (dashes) with E 
in Hartrees. 

depend primarily on the quality of the calculated q = dielectric function, and not on the 
details of the extension to finite q. Using our ab initio dielectric function, we also calculate 
the mean excitation energy and thus stopping powers for relativistic electrons, obtaining 
results in good agreement with experimental data. Furthermore, using the ab initio ADM 
we can extend the stopping power calculation down to energies of O(10) eV, i.e., much lower 
than the Bethe formula, while still maintaining reasonable agreement with experiment. Our 
approach for calculating stopping powers of electrons can be extended to protons and other 
ions by suitably modifying the kinematics. We have also found that the net range (or path 
length) R(E) can be represented as a power law over several decades of energy. That such 
a simple analytic description of the range of probe electrons in a solid exists is interesting 
and also could be useful for reducing computation times in Monte Carlo calculations used 
to study radiation damage and in other applications where numerous electron energy-loss 
tracks need to be considered. In conclusion, we believe our approach has the potential to 
complement or provide an alternative to semi-empirical approaches for calculations of IMFPs 
and stopping powers in condensed matter. 
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